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ABSTRACT 


We  study  homogenization  of  linear  dynamic  thermoelasticity  with 
rapidly  varying  coefficients,  using  a  semi-group  approach.  The 
resulting  homogenized  problem  exhibits  an  unusual  change  in  initial 
temperature . 

A  formal  asymptotic  analysis  predicts  fast  time  oscillations 
in  the  temperature  field.  These  oscillations  explain  the  temperature 
shift,  and  show  that,  for  our  problem,  weak  convergence  in  time  is 
the  best  convergence  that  one  can  obtain. 

INTRODUCTION 

We  discuss  the  problem  of  "homogenizing"  the  equations  of  linear 
thermoelasticity  when  the  mechanical  and  thermal  properties  are  periodic 
and  rapidly  varying.  Following  Bensoussan  -  Lions  -  Papanicolaou  [1] 
and  Sanchez-Palencia  [7]  and  using  a  semi-group  approach,  we  show 
rigorously  that,  as  the  period  of  the  coefficients  goes  to  zero,  the 
solution  of  these  equations  converges  to  the  solution  of  a  related 
constant  coefficient  problem,  the  homogenized  problem.  Then  using 
a  formal  multiple-scales  method,  we  give  what  we  believe  to  be  a 
satisfying  interpretation  of  some  surprizing  features  of  the  results. 

Thermoelastic  behaviour  is  characterized  by  the  coupling  of 
hyperbolic  equations  of  motion  and  a  parabolic  heat  equation.  This 
leads  to  several  interesting  phenomena  in  the  homogenization  process. 

Fast  oscillations  in  the  temperature  field  are  observed;  their  ampli¬ 
tude  remains  finite  as  the  period  goes  to  zero.  Thus  the  solutions 
can  only  converge  in  a  weak  sense  in  time  to  a  slowly  varying  homogenized 


solution. 


Furthermore,  the  initial  data  for  the  homogenized  problem  are 
related  to  the  initial  data  of  the  inhomogeneous  problem  by  a  linear 
transformation  that  is  not  a  projection.  We  know  of  no  other  examples 
of  such  phenomenon. 

In  section  1,  we  formulate  and  prove  the  existence  of  a  homogenized 
thermoelastic  medium.  Section  2  contains  the  more  formal  arguments 
and  the  fast  oscillations  results,  which  are  at  the  root  of  the  ob¬ 
served  change  in  initial  data. 


1.  HOMOGENIZATION  OF  THE  THERMOELASTIC  PROBLEM 

To  cut  down  on  the  overwhelmingly  cumbersome  notations  that 
characterize  thermoelasticity,  we  will  place  ourselves  in  a  scalar 
setting,  that  is  one  where  the  displacement  field  is  taken  to  be 
scalar-valued.  Duvaut-Lions  [2]  show,  using  Korn's  theorem,  that 
this  is  no  loss  of  generality. 

We  consider  a  domain  ft  of  3Rn  .  The  degree  of  smoothness  of 
6ft  will  depend  on  the  type  of  boundary  conditions  adopted.  We  will 
always  assume  that  6ft  is  smooth  enough  for  one  to  be  in  position  to 
apply  Rellich's  theorem  on  compact  imbeddings  of  SoboJev  spaces 
(Folland  [3],  Chapter  6). 

We  will  refer  to  Y  =  II  ]0,  y°t  as  to  the  "reference  cell";  |y| 

i*l 

is  its  volume. 

If  E  is  a  smooth  hypersurface  dividing  Y  into  Y^  and  Y we  define 
aij(y),  Xij(y)  ,  QMy),  £(y)  ,  P (y )  to  be  real  Y-periodic  functions, 
smooth  and  bounded  on  the  closure  of  Y^  and  Y^  but  with  E  as  potential 
surface  of  discontinuity. 
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Furthermore,  ^ (y) ,  A^ ^ (y)  are  assumed 
to  be  symmetric,  strongly  elliptic  on  Y,  that 
is  that  there  exists  a  >  o  such  that: 


(1.1)  a„(y)(resp.  A  _  (y) )  5^  5^  a  ^  on  Y 


3(y)  and  P(y)  are  bounded  away  from  zero.  We  finally  choose  a  such 
that  ot”1  is  a  common  upper  bound  to  the  L^-norms  of  the  coefficients. 

We  extend  all  coefficients  to  all  of  3Rn  by  periodicity.  Our  equations 
are  (Kupradze  [5] ) : 


(1.2) 


~2  e 

pcf»  H- 

£  3t2 


fifx. 
e  (e*  3t 


3 .  ,  3u^~ 

3x.  (aij  (e-  (  3x. 

i  3 


3  t\  3T  x 

3^7  (Aij  (e}  ]  ‘ 

i  J  J 


,x.  £ 

a  (^)  T  )) 
3  b 


a.  .  (*)  a.(-) 
ij  V  jv£; 


3t3x. 

i 


For  the  sake  of  simplicity  we  will  only  consider  Dirichlet  boundary 
conditions  throughout. 


(1-3)  u£  «  0  ,  »  0  on  3Q 


And  for  initial  conditions,  we  will  have: 


(1.4)  uC(x,0)  *  f (x)  ,  (x,0)  =*  g(x)  ,  xe(x,0)  -  k(x) 

e  e 

Our  goal  is  to  study  the  behavior  of  u  and  X  as  £,  the  period, 


2  i  O  «  ;  *  s?  M 
£  1  u  *  I  1  '<♦  5 

i  I  ■  **  '< 


goes  to  zero. 


Vi  v  X  jH 


We  define  H  to  be: 

(1.5)  H  -  x  1^02)  x  L 2(fl) 


On  H,  we  define  the  operator  A^.: 


\ 


(1.6)  A£  = 


3 _  x.  3 _ . 


1  a 4  (?)® 


0.x  ij  e  j  e  3x± 


1  3 _ .  ,x.„  ,x.  . 

P(x)  &xiij  ej  e 


1  4r-<*„&4-> 


/ 


6  A 


with  domain 


(1.7)  D (A_)  =  {u  =  (U,  u„,  X)  CH1^)  x  L,(fl)  x  H^fi) 

t.  t  O  2  .  O 


such  that  A  u  C  H  in  a  weak  sense  } 

t 


Then  the  following  proposition  holds: 


Proposition  1,1 


A£  generates  in  H 


a  strongly  continuous  semi-group  of  operators 


S£(t)  such  that  : 

(1.8)  1 1  s£(t)  1 1  £  a"1  (v  t  >  o  ) 
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Proof 


We  first  consider  for  a  fixed  £  the  norm 


d.9)  |o||  - 


where  ~  denotes  complex  conjugate. 

In  view  of  the  properties  of  the  coefficients,  |  | £  is  indeed 
a  norm  on  H,  equivalent  to  the  natural  Sobolev  norm  on  H,  noted 
||  ||  ,  that  is,  if  0  is  in  H: 


(1.10)  0.  II  oil2  <  |o||<c.-1||u||2 


In  the  norm  |  | £  ,  A£  generates  a  semi-group  of  contractions. 

00 

Indeed,  the  domain  D(A£)  is  dense,  since,  though  Cq  (ft)  functions  do 

oo 

not  belong  to  it,  (ft)  functions  whose  normal  derivatives  are  0  on 
the  only  possible  surfaces  of  discontinuity  for  the  coefficients  (i.e. 
the  ^-scaled  versions  of  £  in  each  of  the  cells  making  up  ft)  do  belong 
to  the  domain  D(A£) .  Checking  that  A£  is  closed,  that  the  range  of 
(1-A£)  is  H  itself  and  that  A£  is  dissipative  offers  no  special  dif¬ 
ficulties  (see  Francfort  [4]  for  full  details).  Note  that  the  measure 
of  the  dissipation, 

(1.11)  Re(A£U,  U)  -  -  Re  </Q  X±j  (t> |^7  dx)  -  ’ 

3  i  2 

(in  view  of  the  properties  of  the  '  s),  is  precisely  the  physical 

dissipation  due  to  heat  fluxes  through  the  boundary. 
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The  result  then  follows  from  the  application  of  Lumer-Phillips's 
theorem  (Yosida  [8],  Chapter  9).  Therefore, 

(1.12)  |se(t)u|e  <_  |u|£  for  any  U  in  H 

and  thus,  using  (1.10), 

(1.13)  ||  Se(t)u||  la"1||u|| 


which  completes  the  proof. 

We  now  leave  the  time  dependent  formulation  and  examine  the 
behavior  of  the  resolvent  of  Ae,  R^  (A£)  as  £  goes  to  0.  At  the  end  of 
this  section  we  will  reintroduce  the  time  dependence  by  using  some  basic 
properties  of  semi-groups. 

It  is  a  direct  consequence  of  (1.8)  (Yosida  [8],  Chapter  9)  that 
the  right  half  complex  plane  belongs  to  the  resolvent  set  of  A£,  for 
every  £.  Let  us  consider  F  =  (f,  g,  k)  to  be  an  element  of  H.  Taking 
A  to  be  real  strictly  positive ,  we  have  the  following  string  of 
equivalences: 


(1.14)  Rx  (A£)F  =  U£  ,  (U£  =  (u 


£  £ 
U*.  ,  I  )  ) 


,  £  £ 

AU  “  ut  *  f 

»<!>  ^  -  V*,t£” 


Ut> 


h-  ‘Vt’lr-1  *  6lt,k  * 

i  j  J  J  i 


P(f)g 
3u* 
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£  £  £ 

The  last  two  equations  (1.15)  have  a  unique  solution  v  «  Xu  ,  T 
1  2 

in  (H  (ft) )  ,  since  the  Dirichlet  form  dP  defined  as: 

o  0 


(1.16) 


de«ve.TE,,(5,m)  -^/saij(f)^-  |i-dx. 


X/n  P<f>vE?  <3*  - 


.x  3C  *  _  /•  >  /X\  3i  3 n 

j  L  dx.  ft  ij  £  3x^  3x. 

i  3  i 


dx 


,x%  e- 


X..,,  /X,  3y 


+  6<t>T  *  d*  +  /q  n 


dx 


1  2 

is  strictly  coercive  on  (HQ(ft))  ,  in  view  of  the  properties  of  the  coef¬ 
ficients. 

£  £ 

If  we  manage  to  find  a  limit  for  u  ,  T  as  £  goes  to  zero,  then 
going  back  up  through  the  string  (1.14)  will  enable  us  to  get  the 
limit  of  R^(A£)F. 

Performing  the  limiting  process  in  (1.14)  is  the  task  of  the 
homogenization  method.  Rather  than  going  through  all  the  lengthy  details 
of  the  argument,  we  merely  mention  the  different  steps  that  were  per¬ 
formed,  underlining  only  the  ones  that  are  not  totally  standard.  For 
further  details  the  reader  is  to  refer  to  Bensoussan-Lions-Papanicolaou 
[1],  Chapter  1*  especially  Sections  3,  9,  and  13,  or,  for  our  problem  to 
Francfort  [4] . 
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tv 


1  2 

First  one  shows  that  u£  and  are  bounded  in  ,  which 

immediately  implies  the  existence  of  a  weakly  convergent  subsequence 
1  2 

in  (Ho (Q) )  converging  to  (u,^).  Since  we  ultimately  show  that  any 
convergent  subsequence  converges  to  the  same  limit  we  do  not  distinguish 
between  the  sequence  and  subsequences  of  this  sequence. 

Then,  defining 


roi  = 


aii (!x -  Uj t£)  the  stress» 

3  j 


(1.17) 


K. 

1 


x  Ail! 

ij  e 


a.  .  (|)a.  A|ii_ 

ID  £  :  e  3xj 


the  heat  flux. 


it  is  easy  to  conclude  that  these  quantities  converge  weakly  in  L2(f2) 
to  (J.  ,  v,  which  in  turn  satisfy: 

f-  2  aai 

pX  u  -  - -  =  P(Xf  +  g) 

I 

(1.18)  < 

_  3*.  _  _ 

*“  a  +  Xv  =  3k  +  a.  ,  ct .  , 

L  3x7  JO  3  3x. 


where  ,  from  now  on,  will  denote  the  Y-average  y^y  dy. 

It  remains  to  determine  CL  ,  and  v.  This  is  the  core  of  homogeni¬ 
zation.  TO  this  effect  we  define  ty) t  'i'(y)  to  be  the  unique 

periodic  solutions,  up  to  a  constant,  in  H^(Y)  of: 
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(1.19) 


3 

WL  (ai3(Y) 


A  °^k 


^  3¥ 

■ST,  lai3ly)  ^ 


S  '  ^Y~  (aij(Y)aj(y)) 


Y  can  be  considered  as  non  standard  with  respect  to  the  "classical 


case .  The  functions : 


(1.20)  wk  =  xk  “  £  X^(e)  '  **  =  *k  '  £  0k(£)  ' 


satisfy: 


x  3w 

^•ij^  ^x7  9x~  dx  =  0 


(1.21) 


3z£  ~~  ^ 

fn  X.  .  (£)  -i  |H_  dx  *  0,  for  any  «*>,U  in  H 
J\l  13  £  3xj 


H^(«) 


Taking  w  and  U  to  be  CMS!)  functions  and  making  use  of  (1.16),  (1.21) 


we  have : 


(1.22)  d£((Xu£,  x£),(ww^,^))  -/Q  a...  If  3^  g- 


x.  ^k  9  «Su  ) 


3zc  ~  £ 

./  x. .(I)-!^1-1^ -  /, op(f)(Xf  +  g) 5 wk dx 

Q  17  £  3xj  3xi  **  £ 
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In  (1.22),  we  have  in  essence  subtracted  from  the  variational 
formulation  of  (1.15)  appropriate  expressions  equal  to  0  in  order  to 
eliminate  products  of  weak  convergences. 

It  is  then  possible  to  go  to  the  limit  in  (1.22)  in  a  way  identical 
to  Bensoussan-Lions-Papanicolaou  [1],  Chapter  1,  Section  3.  Upon 
performing  this  limiting  process  o and  come  out  to  be: 


3u  3X± 

a.  =  (a.  .  -  a.  .  — -  )  Q-  -  (a.  .  a.  -  a.,  a  —i  )T 

x  iD  kD  3y  3x  iD  D  K]  d  3yk 


(1.23) 


=  a  i  39 i  *  3x 

Ki  =  ij  3xk  3x. 


Determining  V  requires  some  extra  effort  and  the  use  of  One 
£ 

defines  g  to  be: 


(1.24)  g£  =  1  +  £  f(f) 


then  it  satisfies,  for  any  ^  in  H  (G) : 


,  x  3g  3u)  ,  ,x.„  .x,  3ui 

(1.25)  !a  (j)  ^7  dx  =  ra  (T>“j <e>  ^  *< 


3  1 


Going  through  the  same  procedure  as  in  (1.22)  but  with  U  equal  to  0 
£  £ 

and  replaced  by  g  ,  we  determine  V  to  be : 


M  ,  3u  ^  .  3Y  .  _ 

(1.26)  V  =  a.  O.  -  a.  -r—  )  - —  +  (a  a,.- —  )  X 

iD  D  3xi  kD  3  3yk 
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o 


Defining  a^ ,  A^,  B^»  X^..,  Cf  to  be 


^  aij  *  aij  *kj  3y 


k 


(1.27)  / 


Ai  *  aijaj  -  \f)  wx 


B.  -  a.  a.  -  ...  |L 
1  il  1  il  3yj 


30. 

X.  .  •  X  .  .  -  X  -r-i- 

il  il  *1  3y 


Y,  -  a.  .a.  -  A. 

i  il  1  i 


a  a.|L 

*1  i 


it  can  be  shown,  using  (1.19) r  that  a^  and  X ^  are  symmetric  positive 
definite,  hence  invertible,  that  and  are  equal  and  that  0  is 
positive. 

We  set: 


(1.28)  Ct±  -  aik  Ak  -  aik  Bk 


Recalling  (1.18),  (1.23),  (1. 26) - (1.28)  yields: 


(1.29)  < 


'*2*  *  *u  -c. 


3f 

.  .  _ - - -  t  a  a .  .«■.  -k —  =  PJC  +  a.  .  iy;-r— 

il  3*.3x.  x]  1  3x.  il  1  3x. 


(F  +  o)Xt  -  X.  .  i— L —  +  Xa.  a.  |H_  =  bk  +  a.  .  (y)a.  (y) 
•  ■  i  3x .  -  - 

i  1 
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( . 


and,  in  view  of  the  properties  of  the  a^'s  and  X^'s,  the  Dirichlet 

1  2 

form  associated  to  (1,29)  is  strictly  coercive  on  (HQ(f2))  ,  hence 

1  o 

(1.29)  admits  a  unique  solution  in  (H  (fl) )  .  Then,  using  (1.14), 

o 

we  end  up  with  the  following  propositions 
Proposition  1.2 


R^(A£)F  converges  weakly  in  (H^(ft))3  to  the  unique  solution  in 
(H1  (fl) )  3  of: 


Xu  -  ut  =  f 

(1.30,  <  Xput  -  a  (|^3—  -W 

1  j  X 


g2x  9ut 

Mia  on  -  ,ijSX*V  5- 

J  1  j  J  J  1 


■5-7—  =  Be  +  Y. 


i  3x. 


We  then  define  A  to  be: 


(1.31)  A  = 


1  ,  3 

-  ij  3x. 3x. 
P  1  j 


&+  a 


d,  .  Ot  ,  11 

i]  :3x7 


-  -  a.  a.  | — 
p  ^  3  9xi 


0  +  a  lj3xi3xj 


It  is  simply  a  matter  of  reproducing  the  proof  of  proposition  1.1, 
but  with  constant  coefficient  this  time  to  show  that  A  generates 
a  semi-group  of  operators  S(t)  such  that: 
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(1.32) 


II  s(t)  ||  <  a' 


for  any  t  >  0 


-1  •  '-1  .  . 
Renaming  ot  the  maximum  of  &  and  &  t  we  deduce  from  proposition 

1.2  and  (1.32)  the  following  corollary. 

Corollary  (1.2) 

1  3 

R^(A£)F  converges  weakly  in  (Hq(&))  to  R^(A)F  where: 


3k  +  Y 


3f 
i  3x. 


(1.33) 


(f,  9'  — rr 


3  +  a 


Now,  (1.8)  implies  that,  for  any  U,  there  is  a  bounded  subsequence 
of  S£(t)U  that  converges  weak-*  in  Ibo(ro+  »  H)  to  G(t)  an  element  of 
I^OR*  ,  H) .  This  is  a  direct  consequence  of  the  separability  of 
L^OR*  ,  H)  and  of  Banach- Alaoglu's  theorem  (Rudin  [8],  Chapter  3). 

Still  identifying  a  sequence  with  its  subsequences ,  we  get  that,  for 
any  V  in  H, 

OO  QO  - 

(1.34)  fQ  e"Xt(S£(t)U,  V)  H  dt  fQ  e'At(G(t),  V)H  dt 


where  (  ,  )H  is  the  natural  inner  product  on  H. 

But  the  resolvent  of  the  generator  of  a  semi-group  applied  on  a  vector 
U  is  equal  to  the  Laplace-transform  of  the  semi-group  acting  on  U 
(Yosida  [8],  Chapter  9)  thus: 


(1.35)  /0°°  e”Xt(Se(t)U,  V)Hdt  -  (Rx(Ae)U,  V)H 
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*  tt 


which  itself  converges  to: 

(1.36)  (RX(A)U,  V)H  =  /0°°  e"Xt(S(t)U,  V)Hdt 


Since  V  is  arbitrary,  we  finally  get,  using  the  uniqueness  of  Laplace 
transforms  of  scalar  function  that: 


(1. 37)  G(t)  -  S(t)U  (t  >  0) 


We  have  proved  in  this  section  the  following  theorem: 

Theorem 

The  generalized  solution  of  (1.2)  with  Dirichlet  boundary  conditions 
and  initial  conditions  (f ,  g,  k)  in  H  converges  weak-*  in  L^UR^  ,  H) 
to  the  generalized  solution  of  : 


(1.38)  < 


3t2 


a.  .  ( 


*2 
3  u 


ij  '  3x.  3x  . 
J  i  1 


-  a 


32t 


V, 


13 


3T  , 

3x. 

x 


-  a. .  a . 
13  3 


a2 
3  u 


3t3x. 


with  Dirichlet  boundary  conditions  and  initial  conditions 


(1.39) 


-  3f 

01  +  Yi  IT 

(f,  g,  — - -  ) 

U  +  a 
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Before  concluding  this  section,  let  us  emphasize  once  more  the 


rather  unusual  change  in  initial  temperature  in  (1.39). 

2.  FAST  OSCILLATIONS  OF  THE  TEMPERATURE  FIELD 

Since,  through  a  weak-*  type  of  convergence  a  rapidly  oscillating 
function  (like  e*^  )  goes  to  0,  it  is  fairly  natural  to  expect  a  ^ 
dependence  of  u  and  x  ,  This  kind  of  problem  is  most  easily  addressed 
using  asymptotic  expansion  techniques.  We  have  already  mentioned  the 
semi-heuristic  character  of  this  section  so  that  we  will  not  dwell 
on  the  restrictions  to  the  problem  that  would  make  the  argument  totally 
rigorous . 

£  £ 

Recalling  (1.2)  we  now  suppose  that  u  and  X  are  functions  of 
k°th  t  and  6=  t/£?  3^  becomes  3^  +  ^  3^.  We  then  Laplace  transform 
(1.2)  with  respect  to  both  t  and  6,  the  dual  variables  being  respectively 
5  and  u  .  From  now  on: 

will  denote  the  t  -  Laplace  -  transform 
^  will  denote  the  6  -  Laplace  -  transform 
v  will  denote  ^  or  ^ 

In  order  to  be  able  to  perform  these  transformations  we  need  to 
inpose  initial  conditions  on  both  t  and  6.  We  will  set: 

u£(x?  0,  5)  *  f(x)  ,  u£(x;  t,  0)  *  p(x,  t) 

e  ^  £ 

“Hx?  0,  6  )  »  g(x)  ,  t,  0)  *  q(x,  t) 

X£(x?  0,  6)  *  k(x)  ,  X£  (x;  t,  0)  *  0(x,  t)  , 
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(2.1) 


where  f,  g,  k  are  as  before  and  p,  q,  0  are  unknown.  We  get: 


P( 


£>{<S2  ~  -  £  )  +  T  (54  Ya£  -  Cp)  + 


4  U 


(2.2)  / 


i  -rt -«)>-}_  <a.j(£><  H-  -  Cj(|)T€  ), 

-  ?>  *  W'c  -  *»>  -  hr  <Vf>!r  > 

^  X  D 


,xx/3  ,rv£  f,  1  3  /\  n 

“  aij(e>aj(e)t^r  (4u  "ir)+e^r^u  -p)} 

J  i  x 


£  £ 

We  seek  an  expansion  of  u  and  t  in  the  form 


(2.3 


i 

Z  £  u^(x,  y,  t,  6) 


*  Z  ^  t,(x,  yr  t,  6)  ,  where  y  =  | 

The  dependence  of  the  u^'s  and  s  on  y  is  taken  to  be  Y-periodic. 
This  is  always  what  is  assumed  when  performing  double  scaling  in  space 
in  problems  related  to  homogenization. 

We  also  need  to  control  the  fast  time  behaviour  of  u^  and  ^ 

Since  we  would  like  them  to  be  oscillating  in  <$,  or,  at  least,  to  be 
such  that 


l  ‘  i  * 

(2.4)  ^  T  -►  +  «-/  u^x,  t,  y,  6)  <j<5  (respectively  T^) 

o 

exist  and  be  finite,  we  are  led,  through  Wiener's  Tauberian  theorem 
(Rudin  16],  Chapter  9)  to  suppose  that: 
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(2.5)  lim  yu-  (respectively  exists  and  is  finite , 


and  we  will  furthermore  assume  that  this  limit  is  to  be  taken  point- 
wise  in  x  and  weakly  in  H1(Y)  with  regard  to  the  y  dependence. 


With  these  considerations  in  mind  we  can  proceed  to  replace 
S  1  3  .  £  e 

by  and  u  and  T  by  their  expansions  in  (2.2) 


3x. 


We  obtain  two  "series”  in  ascending  powers  of  £  starting  at  £ 


-2 


we  identify  the  factors  of  each  of  these  powers  to  0  ,  one  after  another. 
-2 

As  factor  of  £  we  get: 


r 


(2.6)  / 


P(y)  W2u0-Mf-«)  -  ^  (»„(»)  ^  ) 

V 

ax 

(Xij(Y)a^  >  “  aij(Y)0tj(Y)a7"  (uU0  “  =  0 


Since  the  Dirichlet  form  associated  to  the  operator 


(2.7)  D  =  P(y)u2*  -  1^-  (a..(y)  ) 


is  strictly  coercive  on  the  subspace  of  H  (Y)  consisting  of  Y-periodic 

*  1 
functions,  the  first  equation  of  (2.6)  has  a  unique  solution  in  H  (Y) ; 

D  /s 

therefore  ^  ^  is  the  solution.  Hence 

h1  1 . 2 


(2.8>  'uo.£  ♦  a 


But,  in  view  of  (2.5) ,  (2.8)  implies  that  q  *  0  thus  uQ  is  equal  to 

O  V 

S.  and  does  not  depend  on  y.  Inverting  (2.8)  we  get  that  u  does  not 
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o 


depend  on  6  either; 


(2.9)  u  (x,  t)  -  p(x,  t) 


Then  from  the  second  equation  of  (2.6), 


(2.10)  X  -  T  (x,  U)  , 
o  o 


since  the  only  periodic  solution  of  that  equation  is  a  constant  with 
respect  to  y. 

As  factor  of  £  1  we  get,  using  (2.8),  (2.10): 


(2.11) 


Du,  = 


8aij(y^o  3 


1  3Yi  '3x.  'Sy~  ij 


(a. . (y)a.(y)) C 


-  (y)a.(y) 


3  '  3yi 


Defining  ^  to  be  the  unique  periodic  solutions  in  H^(Y)  of 


(2.12) 


2  u  3  xi 

wxi  (akj(y)  ^7“} 


^  -k^wk  1  -  -3 


5yT  1  "  (»kj<y,a3ly,> 


we  obtain  from  the  first  equation  of  (2.11) 


(2.13) 


U  8U0  Hv 

■*! V  * 

3 


Then ,  integrating  the  second  equation  of  (2.11)  with  respect  to  y  and 
defining  and  0^  to  be  the  analogues  of  Y^  and  a  for  and  4^  as  in 
(1.27) , 


18 


‘  o.  . 


+  Y.u  |E- 


(2.14)  T 


x  ax. 

_ _  i 


(def.) 


o  U 


U  +  0^ 


where 


denotes  the  Y-arerage  /  dy 


We  introduce  and  to  be  the  unique  periodic  solutions  in  H^(Y), 
k 


up  to  a  constant ,  of: 


r  3  (\  .  .  v 

-377  (y)  -s—) 


3y^ 


-  “uWjW  37- 


i(zLYP) 
7  k 


(2.15)  / 


^  (Aij {y)  377 } ' 


-  (a.  .(y)ct.(y) 

13  3  *i  I 


The  equations  (2.15)  are  well-posed  since  the  Y-average  of  the  right- 
hand  members  do  vanish  by  definition  of  and  0^  . 

\s 

Recalling  0j (y) ,  we  obtain  for  : 


(2.16)  Ti  =  ”  7  ^3 (y) 3x7  *  Aj 


^  *  •  *  3l^  +  H^1  n  +  an  arbitrary  function  of  x  onlj 


3x . 


Finally,  as  factor  of  £  ,  we  get: 


3u 


rP(Y)  (S2uq  -  Sjj  -  J  +  20*^)  +  Du2  -  |^“  ^aij  ^^"3x7  ]  + 


2 

3  u 


2 

3  u 


3T 


+  a .  .  (v)  - 2.  +  a  (y)- — r — “  -  a.  (y)0t .  (y)  -  -  -3—  (a.  .  (y)°i .  (y)  T  ) 

+  airy;  axi3xj  ijky;3x  3y  ij  y;  Dvy;3xi  By].  ij  3  1 


1  3 


o 


3x, 


ax. 


6(y)  «?„  -  S.  *  u\)  -  |y-  (X. .  <y)^  )  *  §y-<X.  j  (y)^i  ) 


(2.17) 


32t 


32t 


+  Xij  (y)1^7  +  Xij  (y)^I%  "  aij  (y)aj  (y)  3^(C*o  ■  iT  } 

-  ai;j(y)aj(y)  -  a±j  (y)^  (yjlj-cy^)  “  aij  (y)Ctj  (y)^(u^2) 


We  integrate  both  equations  of  (2,17)  with  respect  to  y?  making 
use  of  all  the  previous  results  of  this  section,  we  get: 


(2.18) 


P(C2p  -  Sf  -  g)  +  y3Pu„  * 


2  ij3x^3Xj 


- «  +  *2  *»  ■  ^  %• 


au.  a" 

j-d  a 


3AU  ~2, 


pl.(y)^  IX—  +  UA.  .(y)!2-^ 
ID  3yj  3xi3x]c  ij  3  3yj  aXi 


-c«’f.a.w»-+r7-T‘  1^-^- 

iD  D  3x.  ^(y)U7(y)  3xi 


*i  ,2  tttts 


L  w  aij(y)aj(y)‘5^“  "  u  aij(y>aj(y)^”  ' 


33 i 

3x. 

1 


where  a^j^,  are  to  Xj^  and  ^  what  a^  and  0^  are  to  Xj  ^  in  (1.27)  . 

We  now  consider  the  limit  of  (2.18)  as  U  goes  to  0.  The  following 
result  holds: 

Proposition  2.1 

4^,  U  yHU  go  respectively  to  Xjy  V  #  0  and  0  strongly 
in  H1  (Y) /K  as  y  goes  to  0.  Hence  a  .U  ,<*  M  ,y  w ,0^  go  to  a  ,a  ,Y  ,  a. 

lj  J  1  XJ  J  1 
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o 


The  proof  of  this  proposition,  which  involves  some  basic  esti¬ 
mates  in  H 1(Y)/H  will  not  be  given  here?  refer  to  Francfort  [4] 
for  the  details. 

Proposition  2.1  together  with  (2.5)  is  exactly  what  we  need 
to  perform  the  limiting  process.  Upon  doing  so,  we  come  up  with  a 

A 

set  of  two  equations  for  p  and  H  which,  together  with  the  limit  of 
(2.14),  can  be  interpreted  in  the  time  dependent  domain,  p  and  0 


satisfy: 


P  (3  +  a)n(x,  t)  -  (x,  t) 

0:x,  t)  =* -  * 


(2.19) 


1 


-Bp  Bp  ^  Bn 

P  ^2  *  aij  1^7  -  aijaj  5T 


< 


17  i  cn  a  \ 

18  +  m  Si  xij 

1  D 


-  a.  .a. 


,2 

3_E_ 


ij  j  3t9x. 


Bp 

p(x,  °)  “  f,  3^  ^x'  H(x,  = 


3f 

13k  +  YiBx^ 
3  +  o 


where  0  is  the  limit  of  if  as  |i  goes  to  0. 

It  is  clear  that  r|(x,  t)  can  be  identified  with  I(x,  t)  ,  the  homo 

genized  temperature  field,  and  p(x,  t)  with  u(x,  t) ,  the  homogenized 

displacement  field.  Replacing  0  by  its  value  in  (2.14)  we  also  obtain 

£ 

an  expression  for  the  leading  term  of  the  asymptotic  expantion  of  T  , 
that  is  ^o?  its  6- Laplace  transform  satisfies: 


T 

o 


1  ~g  +  CT 
P  ^ 


n  + 


<V  -  V 


1  +  ^ 


3£_ 

9x. 

1 


(2.20) 
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(2.20)  is  not  invertible  in  general,  but  the  following  proposition  holds 
Proposition  2.2 

and  go  to  zero  as  y  goes  to  +  08 . 

The  proof  of  this  last  proposition  uses  the  same  estimates  as 

the  ones  that  establish  proposition  2.1. 

Propositions  2.1  and  2.2  enable  us  to  conclude  that ,  as  y 

goes  to  0,  T  goes  to  H,  whereas  as  y  goes  to  +  08 ,  goes  to  0. 

In  a  time  dependent  context,  these  facts  translate  into  statements 

on  the  behavior  of  T  near  infinity  and  near  the  origin, 

o 


(2.21) 


<5  -+■  +  00 

< 


^lim 


6-^o  + 


^  ^  Tq(x,  t,  6  )  d<s 

o 

/  T0^x'  t,  6  )  d<S 

o 


t(x,t) 

0  (x,  t) 


The  second  equation  of  (2.21)  is  consistent  with  our  self  imposed 
6  -  initial  conditions.  The  first  equation  shows  that  the  fast 
oscillations  of  the  leading  term  Tq  of  the  asymptotic  expansion  of 

aure  centered  about  T(x,t) ,  the  solution  of  the  homogenized  problem. 


CONCLUSION 

Numerical  evidence  corroborates  the  results  of  section  2  and 
confirms  that  fast  oscillations  are  indeed  the  phenomenon  leading 
to  this  unusual  change  in  initial  data  [4] . 

If  seeking  a  more  physical  explanation ,  one  could  examine  the 
entropy  associated  with  the  problem: 
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*  o 


$(|)T£  +  a 


(*)d  (i)  |y~ 

13  V  jleJ  3x± 


It  is  fairly  straightforward  ,  using  the  results  of  section  2  and 

some  of  the  steps  performed  there,  to  show  that  there  is  no  fast 

dependence  in  time  of  the  space  average  of  the  leading  term  in  the 
e 

expansion  of  s  .  That  the  macroscopic  entropy  of  the  body  is  a 
slowly  varying  quantity  appears  to  be  a  sound  idea  and  does  fit 
our  physical  intuition.  A  fast  oscillation  in  the  temperature  field 
is  the  effect  that  balances  the  space  oscillations  of  the  strains 
due  to  the  inhomogenities  of  the  coefficients  and  allows  the  entropy 
to  evolve  slowly  at  its  own  pace.  In  this  respect  the  unusual 
initial  change  in  temperature  is  needed  to  insure  that  no  fast 
change  in  entropy  is  taking  place  at  time  zero. 

In  contrast  with  other  fast  oscillation  type  problems,  the  "phase" 
of  the  oscillations  in  not  arbitrary  but  perfectly  determined.  It 
also  appears  that  a  geometrical  optics  type  ansatz  in  place  of 
(2.1)  will  fail  since, if  the  solutions  of  (2.12)  are  sums  of  terms 
of  more  than  one  frequency  in  <$,  the  fast  oscillations  need  not  be  periodic 
in  6. 

To  conclude  this  study  let  us  point  out  that  choosing  the 
entropy  as  the  natural  variable  in  place  of  the  temperature  introduces 
space  derivatives  of  the  third  order  and  thereby  prohibits  a  rigorous 
analysis  of  the  type  performed  in  section  1.  A  perturbation  analysis 
using  double  scaling  is  feasible  but  eventually  leads  to  reintro¬ 
ducing  the  temperature  field  as  the  proper  variable. 


23 


ACKNOWLEDGMENT 

This  work  was  supported  in  part  by  ONR  Contract  N00014-76-C-0054 
to  Stanford  University. 


24 


REFERENCES 


[1]  A.  Bensoussan,  J.  L.  Lions  and  G.  Papanicolaou,  Asymptotic  Analysis 
for  Periodic  Structures,  North  Holland  Pub  1978. 

[2]  G.  Duvaut,  J.  L.  Lions,  Variational  Inequalities  in  Mechanics  and 
Physics,  Springer-Verlag  Pub.  1976. 

[3]  G.  Folland,  Introduction  to  Partial  Differential  Equations, 

Princeton  University  Press,  1976. 

[4]  G.  Francfort,  Ph.D.  Dissertation,  Stanford  University,  in  preparation. 

[5]  V.  D.  Kupradze,  Three  Dimensional  Problems  of  the  Mathematical 
Theory  of  Elasticity  and  Thermoelasticity,  North  Holland  Pub., 

1979. 

[6]  W.  Rudin,  Functional  Analysis,  MacGraw-Hill  Publisher,  1973. 

[7]  E.  Sanchez  -  Palencia,  Non  Homogeneous  Media  and  Vibration  Theory, 
Springer-Verlag  Monographs  in  Physics,  1980. 

[d]  K.  Yosida,  Functional  Analysis,  Springer  Verlag  Pub.,  1980,  Sixth  Ed. 


